Variational study of a mobile hole in a two-dimensional quantum antiferromagnet 

using entangled-plaquette states 
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We study the properties of a mobile hole in the t — J model on the square lattice by means of 
variational Monte Carlo simulations based on the entangled-plaquette ansatz. Our energy estimates 
for small lattices reproduce available exact results. We obtain values for the hole energy dispersion 
curve on large lattices in quantitative agreement with earlier findings based on the most reliable 
numerical techniques. Accurate estimates of the hole spectral weight are provided. 

PACS numbers: 02.70.Ss, 71.10.Fd, 75.10.Jm 



I. INTRODUCTION 

Since the fundamental suggestion made by Anderson 
in 1987, the two-dimensional (2D) t - J model (tJM) 
has been regarded as an effective hamiltonian descrip- 
tion of the basic properties of the superconducting cop- 
per oxides. ^'^ Essential physical features of the insulating 
copper-oxide planes of these compounds are successfully 
reproduced by the spin-i antiferromagnetic Heisenberg 
model (AHM),^ to which the tJM reduces when holes are 
absent. On the other hand, such insulating planes can 
turn superconducting if doped with mobile holes, whose 
presence is accounted for in the tJM by adding to the 
AHM a nearest-neighbor hole-hopping term. Hence, the 
t — J Hamiltonian is 



H 



-t > (a 



-h.c.) + J^{S,-Sj-jn,nj) (1) 



where the summations run over nearest-neighbor sites on 
the square lattice, t is the hole hopping amplitude, and 
J the antiferromagnetic coupling between site i and j; 
cfg creates an electron, at site i, of spin projection (e.g., 
along the z axis) s — zL^ avoiding double occupancy. 
Therefore, in terms of the usual fermionic creation oper- 



ator c. 



+ 



,(1 



Si and fii 



with 



Ci,Q, are the spin-i and the number operator 
respectively. 

Despite the simple form of the above Hamiltonian, its 
study constitutes, even in presence of a single hole, a very 
difficult problem. Exact diagonalization (ED), possible 
only for small system sizes (i.e., up to ~ 26 sites), 
approximated analytical treatments, ''"^'^ as well as nu- 
merical approaches of various type,"'^^""'^^ have been largely 
employed to gain insight into the physics of the one-hole 
tJM. Specifically, quantum Monte Carlo (QMC) tech- 
niques based on imaginary time projection furnish, even 
for large systems, accurate (exact in principle) results at 
half filling (i.e., no- hole case) or for t—Q (i.e., static-hole 
case).^'^'^'* However, when t ^ these methods are af- 
fected by the sign problem, for which only essentially un- 
controlled workarounds (e.g., fixed- node approximation) 
exist. Furthermore, being a 2D one, the model of our 



interest cannot be tackled, due to the unfavorable scal- 
ing of the required computational resources, by using the 
Density matrix renormalization group (DMRG),^^ the 
accepted method of choice for ID lattice models. Vari- 
ational Monte Carlo (VMC) schemes do not suffer from 
the sign problem; on the other hand they arc not strictly 
exact, their accuracy depending on the particular choice 
of the trial wave function (WF). For the tJM with a sin- 
gle hole, specific variational WF's^^'^^ have been designed 
affording estimates of the hole dispersion relation in rea- 
sonable agreement with those, generally more accurate, 
obtained by QMC. 

Recently, the suggestion was made of combining 
DMRG and Monte Carlo in order to design novel varia- 
tional algorithms for tensor-network based ansatze ap- 
plicable to systems of spatial dimension higher than 
one.^^'^^ The WF introduced by us,^^ founded on 
entangled-plaquette states (EPS), has been successfully 
employed to study ground-state (GS) properties of lattice 
spin models with open or periodic boundary conditions 
(PBC). Such ansatz has been proven to yield estimates 
of physical observables of accuracy at least comparable 
to that obtainable with the best alternative techniques or 
trial WF's.i®"2i xhc EPS WF can be efficiently optimized 
for plaqucttes of a given size, as well as systematically 
improved by adding more plaquettes and/or increasing 
their size 
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In this paper we apply the EPS ansatz to investigate 
the dynamics of a hole in a quantum antiferromagnet de- 
scribed in Eq. (1). Specifically, we compute the hole en- 
ergy dispersion curve £'(k), as well as the spectral weight 
of the quasi-hole state, for lattice sizes up to iV = 256 
sites in the parameter-range 0.4 < j- < 2.0 (the tJM is 
believed to be relevant to the cuprate superconductors for 
a value of -j ~ 0.4). From the methodological standpoint, 
this study is considerably more challenging than that of 
the various (un) frustrated spin models to which the EPS 
ansatz has been recently applied with success. 

The main goal of this work is to demonstrate the ac- 
curacy of the EPS WF in describing the property of a 
mobile hole in the tJM. This is achieved by a system- 
atic comparison of our results for the energy dispersion 
curve and the hole spectral weight with estimates com- 
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putcd by different approaches. Furthermore, besides the 
mentioned methodological aspects we offer extremely ac- 
curate results for the hole spectral weight in the so called 
"physical region" (i.e.,j ~ 0.4) where the accuracy of 
QMC data, the most reliable ones obtained so far, be- 
comes less satisfactory.^'* Assessing the reliability of the 
EPS ansatz in a case (i.e., the tJM with one mobile hole) 
for which a large number of results obtained with alter- 
native trial states, or other methodologies are available 
in literature, is a fundamental step also in view of fu- 
ture applications of our WF to the long-standing prob- 
lem of understanding the phase diagram of the tJM at 
finite hole concentration. The accuracy of the estimates 
presented here is clearly not a proof of the reliability of 
our variational choice when more holes are present, it 
constitutes, however, the necessary condition to consider 
the EPS WF, with no need of essential adjustments, a 
promising option for the investigation of the finite hole 
concentration scenario. A direct study of the latter prob- 
lem, given the paucity of reliable results obtained in the 
past, would render the assessment of the trustability of 
our ansatz problematic leaving, on the other hand, room 
for dubitative, definitely well-posed, questions such as: 
Why results at finite hole concentration should be taken 
seriously if the reliability of the EPS WF has not even 
been proved in the simpler single hole case ? 

On using a minimal EPS ansatz based on 2 x 2 entan- 
gled plaquettes, we obtain estimates of the hole disper- 
sion curve in qualitative and semi-quantitative agreement 
with those computed by QMC. We achieve quantitative 
agreement, also for other obscrvables, with an EPS WF 
based on plaquettes comprising 9 lattice sites. 

The accuracy of our findings for one hole confirms 
the flexibility of the EPS ansatz, which can be applied 
with no substantial modifications to a large class of 
lattice models. More importantly, the EPS WF used 
in this work is easily generalizable to the case of a fi- 
nite hole concentration, where the sign problem renders 
QMC inapplicable, and the generalization of different 
trial states employed for one or two holes appears quite 
complicated. *^'^^ 

II. VARIATIONAL WAVE FUNCTION 

Let us consider a square lattice of N sites comprising 
one electron per site, a generic WF can be written as 

m=J2win)\n) (2) 

n 

where |n) = n2, . . . , n^), rii is the eigenvalue of erf 
and M^(n) is the weight of a configuration of the system. 
The EPS ansatz is constructed as follows: 

i) Cover the systems with M plaquettes in such a way 
that I sites labelled by ?ii,p,«2,p, ■ • ■ , 't-/,p = np belong 
to the Pth plaquette. 

ii) Express the weight of a global spin configuration as 
a product of variational coefficients C"^ , C"^ , • ■ • , C^f in 



biunivocal correspondence to the particular configuration 
(e.g., spin state along the z axis of the plaquette sites) of 
the plaquettes. Hence, taking into account explicitly the 
bipartite nature of the spin-^ square Heisenberg antifer- 
romagnet at half filling: 

(n|$) = W[n) = (-l)^(")7'(n) (3) 

being F(n), according to the Marshall-sign rule, the num- 
ber of "down" spins in one of the two sublattices^^ and 

V{n) = n Cr (4) 

p=i 

the crucial part of the ansatz. 

Analogously, in presence of a single hole, our EPS WF 
is defined via 

(mlV'k) = W{m) = e-*"^^" {-lf^"'^P{m); (5) 

here m refers to the generic configuration containing one 
hole and A^ — 1 electrons, k and Rh are the hole mo- 
mentum and the spatial position of the empty site re- 
spectively. It is worth mentioning that the trial states 
adopted in this work naturally incorporate, for plaquettes 
including more than one site, electron-electron as well as 
electron-hole correlations. We estimate GS energies of 
the Hamiltonian in Eq. (1) via variational optimization 
of the states (3) at half filling and (5) in presence of a 
single hole. Independent optimizations are carried on for 
each k value by Monte Carlo sampling both the energy 
and its derivatives with respect to the varational param- 
eters (i.e., the Cp'"s). Further details concerning the 
variational family of EPS and the numerical technique 
adopted to minimize the energy are given in Refs. 18 
and 20. Although the EPS ansatz is exact when a single 
plaquette as large as the system is employed, accurate 
results are already obtainable by using N entangled (i.e., 
overlapping) plaquettes comprising 4 sites, the accuracy 
of the estimates being sensibly improvable increasing the 
plaquette size. Specifically, we compute the hole energy 
dispersion curve 

E{k) = E"^{k)- E^^"^, (6) 

defined by the difference between the k-dependcnt GS 
energy of the one-hole tJM and that of the AHM, and 
the hole spectral weight 

z(k) = |(V^k|ck|<i>)r, (7) 

where jV'k) (1*1')) is the normalized GS of the one (zero) 
-hole system. 

We compare the EPS results with exact ones for a 4 x 4 
lattice and with the most accurate available in literature, 
for larger systems. In this case the main features of the 
hole dispersion relation, namely the characteristic shape 
with a global minimum at k — (±^,±^) and a nearly 
fiat region around k = (±7r,0), are very well reproduced 
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with an ansatz based on 2 x 2 plaquettes, while the in- 
creasing of the plaquette size to 3 x 3 gives rise to an 
approximately rigid shift, towards lower energies, of the 
band and an almost perfect agreement (not reachable, to 
our knowledge, with other variational WF's) with previ- 
ous QMC studies for all the j values considered in this 
work. Moreover, with the 3x3 EPS, accurate estimates 
of the hole spectral weight arc obtained, on the same 
footing; via the Monte Carlo method. 



III. RESULTS 

The GS energies of the one-hole tJM on a 4 x 4 lattice 
can be obtained, given the small system size, by using a 
single plaquette comprising all the lattice sites. In this 
case the EPS ansatz is exact, the error bar of our esti- 
mates being only due to the limited number of configura- 
tions sampled, and reproduces, as expected, ED results.^ 

Clearly, for larger systems, the use of a single plaque- 
tte is not feasible, therefore, it is important to asses the 
accuracy of our WF for smaller plaquette size. The error 
relative to the exact result: E"^{±§,±§) = -21.161 
of our estimate of the global energy minimum at -j = 0.5 
is ~ 3% for the 2x2 EPS ansatz, decreasing down to 
less than 0.2% when 3x3 plaquettes are used. For this 
plaquette size, our results compare favorably to those 
obtained with different trial states, even when larger 
lattices are considered. For example, the EPS energy 
when k = (±|,±|) on a 8 x 8 lattice at j = 0.4 is 
—78.463(7) J, while the WF including both spin-spin and 
momentum dependent spin-hole correlations adopted in 
Ref. 11 yields a value of -76.753(14)J (i.e., ~ 2% 
higher) . 

As opposed to the case in which more holes are present 




(0,0) (7t,0) (7t,7t) 

k 



(0,0) 



FIG. 2: (color online). Hole energy dispersion curve for a 
lattice of 256 sites with PBC at j = 0.8. EPS estimates 
(triangles) are computed using square plaquettes of 9 sites. 
QMC results (crosses) are also shown for comparison.^* 



and the sign problem is particularly severe, QMC tech- 
niques based on imaginary time projection, for the one- 
hole tJM, provide extremely accurate results, to our 
knowledge the most accurate in literature, when, due to 
the large system size, the ED of the Hamiltonian is not 
possible. Hence, a comparison of the EPS WF against 
these approaches is in order. 

The hole energy dispersion curve, defined in Eq. (6), 
is shown in Fig. 1 for the 16 x 16 lattice at j = 2.0. 
Our estimates (triangles), computed with the 3x3 EPS 
ansatz, are indistinguishable (taking into account the er- 
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FIG. 1: (color online). Hole energy dispersion curve for a 
lattice of 256 sites with PBC at f = 2.0. EPS estimates 
(triangles) are computed using square plaquettes of 9 sites. 
QMC results (crosses) are also shown for comparison.^* 



FIG. 3: (color online). Hole energy dispersion curve for a 
lattice of 256 sites with PBC at ^ = 0.4. EPS estimates 
are computed using square plaquettes of 4 (circles) or 9 (tri- 
angles) sites. GFMC results (crosses) ^"^ are also shown for 
comparison. 
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TABLE I: Band width (W) and energy difference (G) between 
E{n,0) and E{^, ^) , in units of J, for a system of 256 sites 
with PBC. Error bars are in parenthesis. 



J 

t. 


W 


G 


2.0 


0.53(3) 


0.15(3) 


0.8 


1.54(3) 


0.39(3) 


0.4 


2.63(4) 


0.74(8) 



ror bars) from available QMC oncs^^ (crosses). Around 
k = (tTjO) the curve is nearly flat, in agreement with 
other studies. ^^'^'^ It has to be mentioned that the de- 
generacy observed on the 4x4 lattice between E^tt, 0) 
and £'(f,|) is due to geometrical reasons^ and is no 
longer present when the lattice size increases. 

The hole dispersion relation for lower values of j is 
plotted in Figs. 2 and 3. The main qualitative fea- 
tures of the curve (e.g., the position of the minimum) 
appear independent of the value of j . To illustrate fur- 
ther the accuracy of the EPS WF, estimates of i?(k)/t 
obtained with 2x2 plaquettes (circles in Fig. 3) are also 
shown for j = 0.4. Despite the relatively small num- 
ber of variational coefficients utilized for this ansatz, the 
resulting dispersion relation is in close agreement with 
that found with 3x3 plaquettes as well as with Green 
function Monte Carlo (GFMC) results^"^ (triangles and 
crosses respectively in the same figure) and QMC ones.^^ 
Specifically, the hole band predicted by using such a min- 
imal EPS ansatz differs from the more accurate predic- 
tion achievable with 3x3 plaquettes only by a rigid en- 
ergy shift towards lower energy values. This energy shift 
(albeit quantitatively different) occurs for all the three 
values of j considered here, pointing out how the EPS 
ansatz, even in its simplest and economical implementa- 
tion (i.e., when 2x2 plaquettes are used), can provide 
reliable estimates of observable such as the effective mass, 
the band width, and generally all those obtainable from 
band energy differences. 

Values of the band width W and the energy gap G 
respectively defined as = £^(0,0) - E{^, f ) and G = 
E{tt^ 0) — £'(§, §) are reported, as a function of j, for a 
16 X 16 lattice in Tab. I. Our results of the band width 
are m excellent agreement with QMC^** and GFMC" 
calculations. At j = 2.0 our estimate differs from the 
VMC one: 2.11(3) reported in Ref. 12. Moreover we 
find that G, at least in the cases examined here, varies 



TABLE II: Hole spectral weight at k = (|, |) for a 16 x 16 
lattice with PBC. Estimates obtained with VMC based on 
a different WF,^^ QMC," and SCBA^" are also shown for 
comparison. Error bars are in parenthesis. 



J 
t 


EPS 


VMC 


QMC 


SCBA 


2.0 


0.596(2) 


0.663(3) 


0.58(4) 




0.8 


0.427(3) 




0.40(4) 


0.504 


0.4 


0.340(2) 


0.375(2) 


0.32(2) 


0.34 
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FIG. 4: (color online). Hole spectral weight for a 16 x 16 
lattice with PBC at j — 0.4. The dashed line is a guide to 
the eye. 



linearly with t according to the relation j- ^ 0.3. 

Next we discuss our findings for the hole spectral 
weight Z (sec Eq. (7)). On a 4 x 4 lattice wc obtain, for 
{ = 0.4, Z(f , f ) = 0.3996(5), in agreement with the ED 
value'* of 0.4. EPS values of Z(|, |) as a function of j- 
on a 16 X 16 lattice are shown in Tab. II. Our estimates, 
computed with plaquettes of 9 sites, are in agreement 
with QMC ones^^ for all the values of j; agreement is 
found with self consistent Born approximation (SCBA) 
calculations"'^° at j = 0.4, this approach providing, as 
well as the WF used in Ref. 12, higher values of the 
spectral weight when j increases. 

The momentum dependence of the spectral weight for 
J = 0.4 is shown in Fig. 4 along particular cuts of the 
Brillouin zone. ZCk) is approximately 13% at k = (0, 0) 
and, going in the (1,0) direction, reaches ~ 0.4 at k = 
(tt, 0); it then decreases, slightly, down to ~ 0.34 at k = 
(■|,^), and more pronouncedly, in the (1,1) direction, 
beyond this k value. 



IV. CONCLUSIONS AND OUTLOOK 

The two dimensional t — J model in presence of a sin- 
gle mobile hole has been investigated by means of vari- 
ational Monte Carlo simulations based on the EPS trial 
wave function. For systems comprising up to V = 256 
sites (i.e., the maximum size considered in this work), an 
EPS ansatz consisting of N entangled plaquettes of 4 sites 
provides estimates of the dispersion curve in qualitative 
and semi-quantitative agreement (up to a rigid energy 
shift) with the most accurate alternative numerical ap- 
proaches. By increasing the plaquette size to 9 sites, with 
a minimal additional computational effort, quantitative 
agreement is easily recovered and accurate results for the 
hole spectral weight obtained. 
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Given the accuracy of the results presented in this 
work, as weh as the simphcity of the variational algo- 
rithm used to minimize the energy, the EPS ansatz ap- 
pears a very promising option for future investigations of 
the t — J model in the case of finite hole concentration, a 
problem of great physical interest to which exact quan- 
tum Monte Carlo techniques based on imaginary time 
projection are not applicable due to the sign problem. 
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